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We suggest a practical way to estimate the condensate fraction of an interacting dilute Bose gas 
confined by an external harmonic potential as a function of temperature and scattering length. 
It shows that an increase of the scattering length produces an exponential decrease of condensate 
fraction. 
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Over the last decade the phenomenon of Bose-Einstein 
condensation of ultra-cold dilute alkali Bose gas confined 
in a magneto-optic trap has been a subject of fundamen- 
tal interest. Even below the transition temperature some 
atoms leave the lowest energy state and remain in the ex- 
cited state by an atom-atom interaction. This causes a 
depletion of the condensate number in comparison to its 
ground state occupation for the non-interacting case, and 
it depends on temperature and interaction strength both. 

The study of the condensate fraction (CF), Nq/N, 
where N being total number of atoms, of a dilute Bose 
system has a long history. The first microscopic theory 
was suggested by Bogoliubov. If the system is homoge- 
neous and dilute (pa 3 <C 1) where p is the density, the 
CF at T = OK is known from the Bogoliubov theory as 
1 — P fl3 d 2]. However, this formula is not appli- 

cable for an inhomogeneous system because confinement 
makes the wave function of the Bosons space- varying and 
so is the density inhomogeneous. 

In this paper we suggest a practical way to estimate 
the CF of interacting Bosons in a harmonic trap. Our 
approach is to start with a non-interacting harmonically 
trapped Bosons at T ^ and to extend the argument for 
the interacting case. Some experimental variables will be 
applied for the practical purpose. 

Consider a three-dimensional(3D) system of non- 
interacting bosons confined by a harmonic potential of 
angular frequency uj. Then, the number distribution 
function in the grand canonical ensemble at temperature 
T is given by 
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where (3 — 1/fcgT and p is the chemical potential. e n 
is the discrete energy spectrum of the harmonic oscilla- 
tor given by e n — h(n x u> x + n y u y + n z uj z ) + eq, where 



0, 1, 2, ... and £q is the zero point energy. 



We rewrite Eq. (fTJ as a sum of two parts: ground state 



and excited states as 
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where d(n) = (n + l)(n + 2)/2 is the degeneracy of the 
harmonic oscillator. If we consider an isotropic trap of 
lu x = uj y = uj z = uj, the number of atoms in the conden- 
sate at temperature T is well-known as [3, 0, Q 
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where z is the effective fugacity defined by z 
and g s is the Bose series function defined by g s (z) — 
^^L 1 z n /n s . Eq. (|3|) is solved numerically to obtain 
N^ on (T). Note that the transition temperature for N = 
10 4 trapped atoms with uj = 10 3 sec~ 1 is about 150nK. 

In the interacting system of the interaction energy in 
the ground state U [No, a), the zero point energy is shifted 
from Sq to Eq + U . Then, it depends on number of con- 
densate atoms No and the s-wave scattering length a, 
both. In this scheme, the new ground state occupation 
number can be rewritten as 
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The ground state information of a dilute Bose gas is 
obtained from the solution of Gross-Pitaevskii equation 
as 
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where Uq(= 47rft 2 a/m) is the interaction strength be- 
tween atoms and V^o is the macroscopic wave function 
of the condensate which is normalized to the number of 
particles in the ground state, No- The macroscopic wave 
function *o is obtained from the solution of Eq. ([5]) nu- 
merically and then, the ground state interaction energy 
U(N ,a) is given by [1, 0| 
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FIG. 1: Interaction energy per atom of harmonically trapped 
and interacting Bose atoms. The unit is flu. 
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FIG. 2: Condensate number Nq as a function of scattering 
length shown for various temperatures. 



We plotted the interaction energy in Fig. 2 as a func- 
tion of a dimensionless interaction strength g(No,a) — 
A^N^a/ 'a/ lo , where a^o is the harmonic oscillator length 



given by a^o — \fhJmuo. The U is not linear but close to 
a logarithmic function for small g [![. 

A direct calculation of Eq. is not reliable because 
Eq. ([3]) is too sensitive to get an accurate z. Instead, 
comparing Eq. (j4]) with Eq. ([2]), we may obtain a prac- 
tical form of the CF as 
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This approximation is valid only when \(3U\ < 1. In the 
non-interacting limit Nq — > Nq 071 , and otherwise No falls 
exponentially as U increases. 

Because U itself is a function of No, the CF has to be 
obtained self-consistently. We obtain Nq oti (T) from Eq. 
(J3j) and then put it into g(N ,a) to solve the Eq. (|5|) 
numerically. Then, substituting the solution of Eq. (|5|) 
into Eq. ^ again, we obtain the interaction energy U 
as a function of No and a. Next, substitution of U into 
Eq. (|7|) again, we obtain a new iVo- Then, it is put back 
into Eq. ([5|) for self-consistency. 

Calculations are carried out for \/3U\ < 0.1 with the 
practical variables of N = 10 4 and ui = 10 3 sec~ 1 . The 
results are plotted in Fig. 3. The Nq shows an exponen- 
tial decrease as scattering length a increases. The slope 
is steeper for lower temperatures. As scattering length 
increases with other parameters remaining fixed, the de- 
pletion of the condensate is enhanced. 

We suggested a practical way to estimate the conden- 
sate fraction of interacting and trapped Bose atoms in 
specific experimental ranges. The effect of the attrac- 
tive interaction is not clear because it may increase CF 
a little bit for negative scattering length compared with 
non-interacting system. That is, there exists a possibil- 
ity an attractive interaction between fermions and bosons 
may increase CF compared with that of bosons only. 

The author thanks Mukunda P. Das and G. Gnanapra- 
gasam for useful discussions. 
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